Insulating phases of electrons on a zigzag strip in the orbital magnetic field 
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We consider electrons on a two-leg triangular ladder at half-filling and in an orbital magnetic field. In a two- 
band regime in the absence of the field, the electronic system remains conducting for weak interactions since 
there is no four-fermion Umklapp term. We find that in the presence of the orbital field there is a four-fermion 
Umklapp and it is always relevant for repulsive interactions. Thus in this special ladder, the combination of the 
orbital magnetic field and interactions provides a mechanism to drive metal-insulator transition already at weak 
coupling. We discuss properties of the possible resulting phases C0S2 and various COSl and COSO. 

PACS numbers: 



I. INTRODUCTION 

This paper complements our earlier work Ref . [l| on the ef- 
fects of Zeeman field on a Spin Bose-Metal (SBM) phased 
(the reader is referred to Refs. \V^ for general introduction). 
Here we consider the orbital magnetic field on the electronic 
two-leg triangular ladder 

Previous studies of ladders with orbital field were done on a 
square 2-leg case and mainly focused on generic density (see 
Refs. [j-la and citations therein), while the triangular 2-leg case 
has not been considered so far. In the context of Mott in- 
sulators at half-filling, microscopic orbital fields were shown 
to give rise to interesting scalar chirality terms operating on 
triangles in the effective spin Hamiltonianii— On the other 
hand, it was also arguedii— that if a Mott insulator develops 
a noncoplanar magnetic order with nontrivial chiralities, this 
can imply spontaneous orbital electronic currents. 

In this paper, we focus on the simplest ladder model with 
triangles, the zigzag strip, and discuss instabilities due to ex- 
istence of orbital magnetic field and properties of the result- 
ing phases. Our main findings are presented as follows. In 
Sec. ini we determine the electron dispersion in the orbital 
field and perform weak coupling renormalization group (RG) 
analysis in a two-band regimcia^— Unlike the case with no 
field, we find that there is a four-fermion Umklapp interaction 
which is always relevant for repulsively interacting electrons 
and provides a mechanism to drive the metal-insulator transi- 
tion. This Umklapp gaps out all charge modes and produces 
a C0S2 state. In Sec.|lll]we describe physical observables in 
this phase, and in Sec.|IV]we analyze possible further instabil- 
ities in the spin sector and properties of the resulting phases. 
We conclude with discussion of the orbital field effects in the 
context of the Spin Bose-Metal phase of Ref. U where the Mott 
insulator is first produced by an eight-fermion Umklapp and 
the new four-fermion Umklapp appears as a residual interac- 
tion. 
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FIG. 1: Top: Zigzag strip with uniform flux $ penetrating each tri- 
angular plaquette. Bottom: Convenient representation of the model 
as a ID chain with first- and second-neighbor hoppings. We choose 
a gauge such that tx,x+i = ii and tx,x+2 = ^26'*''°'' '""""'. The unit 
cell consists of two sites labelled A and B. 



the triangular strip with uniform flux <I> passing through each 
triangle. Figure [T] illustrates our gauge choice. 
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Here and throughout, we refer to sites by their ID chain co- 
ordinate X. Since the second-neighbor hopping depends on 
whether x is even or odd, the unit cell has two sites which 
we label A and B. The Hamiltonian for such an interacting 
electron system is H = Hq + Hy, with 



Hn = - 
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II. WEAK COUPLING APPROACH TO ELECTRONS ON 
A ZIGZAG STRIP WITH ORBITAL FIELD 

Let us apply weak coupling Renormalization Group (RG) 
to study effects of electronic interactions in the presence of the 
orbital magnetic field. We start with free electrons hopping on 



In the first and last lines, we suppressed the sublattice labels, 
and ?i(x) = J2a '^ai^) ^01 (x)- We assume that Hy is small 
and treat it as a perturbation to Hq. The free electron disper- 
sion is 

e(fc) = ±2y/[ti cos (fc)]2 + [t2 sin ($) sin {2k)]^ 

-2*2 cos ($) cos (2fc) - ^I . (7) 




FIG. 2: Free electron spectrum in the presence of the orbital field, 
cf. Fig. [T] Here ^(fc) is given by Eq. Q with two branches and we 
focus on the regime when both bands are partially populated; we 
take ti = 1, i2 = 1, and $ = vr/lOO for illustration. The half-filling 
condition requires kpi + kF2 = 7r/2. 



We are focusing on the regime with two partially filled bands 
as shown in Fig. |2] For small flux, this regime appears 
when t2/ti > 0.5. We denote Fermi wavevectors for the 
right-moving electrons as kpi and kp2 and the correspond- 
ing Fermi velocities as fi and V2- The half-filling condition 
reads kpi + kp2 = tt/2. 

The electron operators are expanded in terms of continuum 
fields. 
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where P = R/L = +/— denotes the right and left movers, 
a = 1,2 denotes the two Fermi seas, and M = A or B de- 
notes the sublattices. In the specific gauge, the wavefunctions 
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with 

{sin(6'fe), cos(6'fc)} ex {ii cos(fc), is sin($) sin(2fc)}(10) 

Note that k belongs to the reduced Brillouin zone 
[-7T/2,7r/2]. 

Few words about physical symmetries. The present prob- 
lem has SU(2) spin rotation symmetry (7?.) but lacks time re- 
versal because of the orbital field. It also lacks inversion sym- 
metry and translation by one lattice spacing. However, the 
system is invariant under combined transformations such as 
inversion plus complex conjugation (/* : x — > — x, i -^ — i) 
and translation by one lattice spacing plus complex conjuga- 
tion (Tj* : a; — )• a; + 1, i -^ ~i). Table Ulists transformation 



TABLE I: Transformation properties of the continuum fields under 
I* (inversion plus complex conjugation), Ti (translation by one lat- 
tice spacing plus complex conjugation), and TZ (SU(2) spin rotation 
about arbitrary axis n by an angle (f>). We also show transformation 
properties of bilinears E\.2 defined in Eqs. J4U-ll42t. 
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properties of the continuum fields under these two discrete 
transformations and under the SU(2) spin rotation. Since the 
symmetries are reduced compared to the case without the or- 
bital fieldr^^ we need to scrutinize interactions allowed in 
the continuum field theory. 

Using symmetry considerations, we can write down the 
general form of the four-fermion interactions which mix the 
right and left moving fields: 
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where we defined 
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Note that besides the familiar momentum-conserving four- 
fermion interactions H^ and 'H'^, there is also an Umklapp- 
type interaction "H". 

Using the symmetries of the problem, we can check that all 
couplings are real and satisfy wi2 = W21 and A12 = A21, and 
we also use convention wu = W22 = 0. Thus there are 9 
independent couplings: W12 ' •^ii'^i -^22 ' -^12 ' ^'i'^ "4- 

With all terms defined above, we can derive weak coupling 
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Here O = dO/di, where £ is logarithm of the length scale. 
We have also defined 
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We can obtain bare values of the couplings for any elec- 
tronic interactions by expanding in terms of the continuum 
fields. In the case of small flux, the couplings \pI" and w''l'^ 
in Eqs. (fTTl i and (fT2] i are only modified slightly and can be 
treated as the same as in Ref. [l^ with extended repulsion. 
For the coupling u^ in Eq. ([Til l, the bare value of u^ in the 
small flux limit is Yl,x' ^{^ " x')e't(^-^') x |^[sin(fci?i) + 
sin(/c_F2)]<l' oc $, where x and x' belong to the same sublattice 
(A or B). Therefore, we can see that the parameter U4 which 
measures the strength of the umklapp process is linearly pro- 
portional to the flux and goes to zero if we gradually switch 
off the flux. For repulsive interactions, we generally expect 
positive A'' (see, e.g., Ref. T4 with extended repulsion). Then 
according to the RG Eq. (|24] |. positive initial Aj2 will drive u^ 
to increase exponentially. Thus we conclude that the starting 
two-band metallic phase is unstable due to the new Umklapp 
term. 

To analyze the resulting phase(s), we use bosonization to 
rewrite fermionic fields in terms of bosonic fields, 
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with canonically conjugate boson fields: 
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ITrSabSaP Q(x~ x) , (28) 



where Q{x) is the Heaviside step function. Here we use Majo- 
rana fermions {riaayVbp} = 26abSap as Klein factors, which 



assure that the fermion fields with different flavors anticom- 
mute. 

It is convenient to introduce new variables 
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(31) 



and similarly for ip variables. We can then write compactly all 
nonlinear potentials obtained upon bosonization of the four- 
fermion interactions: 

H" = 4u4f sin(2¥>p_) sin(26'p+) , (32) 

W = {w1^Jri2Jli2 ~ wIJri2 ■ JL12) + H.c. = (33) 
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m^vnm^mi 



(36) 
(37) 

(38) 
(39) 



We will not analyze the RG flows in all cases. Our main 
interest is in exploring the orbital magnetic field effects on 
the C2S2 metallic phase and nearby Cl[p— ]S2 spin liquid. 
Therefore we consider the situation where in the absence of 
the W4 term we have the stable C2S2 phase described by RG 
flows such that A^^ reach some fixed point values, w'^i^" are 
irrelevant, and A^^ are marginally irrelevant - this is realized, 
for example, in Ref. [13 for sufficiently long -ranged repulsion. 

As we have already discussed, for repulsive interactions 
we expect X'^^ > and hence any non-zero U4 will increase 
quickly. In this setting it is then natural to focus on the effects 
of the H"" first. From the bosonized form Eq. ( l32b . we see that 
it pins 



sin(20p-i_) =: — sign(u4) sin(2(^p_) = ±1 



(40) 



Thus, both "p^" and "p+" modes become gapped and the 
system is an insulator. This insulator arises because of the 
combined localizing effects of the orbital field and repulsive 
interactions. 

Having concluded that U4 becomes large, if we were to con- 
tinue using the weak coupling RG Eqs. (fT6l)-(l24b. we would 
find that U4 drives A12 to large positive value, which in turn 



drives A'' to negative values and destabilizes couplings u'12 , 
and all couplings eventually diverge. If we do not make finer 
distinctions as to which couplings diverge faster, we would 
conclude that the ultimate outcome is a fully gapped COSO. 
We will analyze different COSO phases arising from the com- 
bined effects of M4 and A"^ later Here we only note that the 
bosonized theory suggests that a C0S2 phase can in principle 
be stable. Indeed, once we pin tpp_ to satisfy Eq. (|40] |. the W 
interaction vanishes leaving only the effective A*^ couplings 
in the spin sector. The stability in the spin sector is then de- 
termined by the signs of the A"^ couplings. If A^^ > 0, the 
spin sector is stable and we have the C0S2 phase. In what 
follows, we will identify all interesting physical observables 
in this phase and will use it as a starting point for analysis 
of possible further instabilities and features of the resulting 
phases. 



III. OBSERVABLES IN THE MOTT-INSULATING PHASE 
IN ORBITAL FIELD 

To characterize the induced insulating phase(s), we con- 
sider observables constructed out of the fermion fields. The 
only important bilinear operators are 
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and their Hermitian conjugates. All other bilinears contain 
field 9p- and hence have exponentially decaying correlations 
once (fip- is pinned. Here and below, repeated spin indices 
imply summation. Operators Ei , E2 are scalars and Vi , V2 
are vectors under spin SU(2). One can check that Ei and E2 
have identical transformation properties under all symmetries 
and therefore are not independent observables, and the same 
holds for Vi and V2 ■ 

The scalar bilinears Ei and E2 appear, e.g., when express- 
ing fermion hopping energies and currents. Specifically, con- 
sider a bond [x, x' = x + n] (we will focus on n = 1 or 2), 

i3("'(a;) - t^^^+nciix)caix + n) + H.c. , 



J^(")(x) - i[t^,^^^+ricl{x)ca{x + n) - H.C.; 



(45) 
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where we have suppressed "sublattice" site labels A or B and 
tx,x+n is defined in Eqs. ([T]!-©. In general, we need to con- 
sider separately cases [x G A,x' £ A], [x £ B,x' G B], 
[x G A,x' £ B], [x G B, x' S A]. After expansion in terms 
of the continuum fields in each case, we find that all cases can 
be summarized by a single form that requires only the physi- 
cal coordinate x but not the sublattice labels: 
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where A 



1,2,3,4 



are some real numbers. The above concise 



form is possible because of the T^ symmetry involving trans- 
lation by one lattice spacing. 

In our analysis below, we will also use a scalar spin chirality 
defined as 



Xix) = S{x) ■ [Six - I) X Six + I)] 



(49) 



From the perspective of symmetry transformation properties, 
the scalar spin chirality and the so-called 'site-centered' cur- 
rents 



xix), J^^\x-1), J^'\x - I) + J^'\x) 



(50) 



have the same transformation properties. (Note that the above 
currents are named site-centered because they get inverted un- 
der inversion about site x. Similarly, we can also call J''-^^x) 
to be 'bond-centered' since it is inverted under inversion about 
X + 1/2, the center of the bond between x and .t + 1.) 

Thus, up to some real factors, we can deduce that the scalar 
spin chirality in Eq. ( |49] l contains the following contributions 
(focusing on terms that have power law correlations): 
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The vector bilinears Vi and V2 appear when expressing spin 
operator. 
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We consider separately two cases x E A and x G B. After 
expanding in terms of the continuum fields, we find that both 
cases can be summarized by a single form that requires only 
the physical coordinate x. 
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Six) 



where A[ 2 are some real factors. 

The bosonized expressions for E12 are: 
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The bosonized expressions for Vi and V2 are similarly 
straightforward. Since we have SU(2) spin invariance, for 
simplicity, we only write out V^ : 
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Besides the bilinears considered above, we have also iden- 
tified important four-fermion operators, 



B. 



(1) 

stagg 



J ^0 



t „o. 



I = *(Cm'^ Cli)(c)j2^ Cl2)+H.C. 



cos(20<x+) +cos(26lff_) sin(26'p+) 



(58) 



B. 



(1) 

stagg 



II = i(c|,^iCTCLl)- (c}j2i?CL2)+H.C. 



^. 



stagg, I 



qZ 

•^staggJI 



cos(26i^+) - cos(26i<,_) + 2rcos(2(p<,_) x 
X sm{20p+) ; (59) 
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a above is the 2x2 identity matrix and a are the usual Pauli 
matrices. The label "staggered" informs how they contribute 
to the spin and bond energy observables. 
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As an example, the above contributions to the bond energy 
arise from expanding nearest-neighbor energies n{x)n{x + 1) 
and S{x) ■ S{x + 1) in terms of the continuum fields. Again, 
we need to consider separately cases x G A or x € B, but 
we find that both can be summarized by the form that requires 
only the physical coordinate x. 

Note that we have only listed observables containing 
sin (26'p_|-). Expressions that contain cos(20p_|-) vanish be- 
cause of the pinning condition Eq. ( |40l ); in particular, there 
is no B^~J'^ . Also, for brevity we have only listed the 
bosonized form of the z-component of the spin observable. 

There are several other non-vanishing four-fermion terms. 
Thus, there is a term which can be interpreted as a staggered 
scalar spin chirality; however, it is identical to H", Eq. ( fTSl ). 
and is always present as a static background in our system. In 
addition, there is a spin-1 observable which can be interpreted 
as a spin current, and a spin-2 (i.e., spin-nematic) observable. 
In the C0S2 phase, these will have the same power laws as 
'^stagg ^iid "S'stagg- However, in our model, they become short- 



ranged if any spin mode gets gapped, and we do not list them 
explicitly as the main observables. 

Let us briefly describe treatment of the Klein factors (see, 
e.g., Ref. [13 for more details). We need this in the next sec- 
tion when determining "order parameters" of various phases 
obtained as instabilities of the C0S2 phase. The operator 
r = ?7it?7i4,?/2t'?2i has eigenvalues ±1. For concreteness, we 
work with the eigenstate corresponding to +1: r|+) = |+). 
We then find the following relation 
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TABLE II: Summary of the properties of the phases from different 
instabilities in the spin sector. 



and the scalar bilinears are expressed as 
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For repulsively interacting electrons, the Umklapp term "H" 
appearing in the presence of the orbital field is always relevant 
and pins 9p^ and (pp_ as in Eq. ( |40l i. As already discussed, 
for such pinning the VF-term Eq. ( |33] | vanishes. Therefore, 
as far as further instabilities of this C0S2 Mott insulator are 
concerned, we need to discuss the Vj_ -terms Eq. (|38] | that can 
gap out fields in the spin sector. 

The instabilities depend on the signs of the couplings A^j, 
A22, and Aj2, so there are eight cases. The simplest case is 
when all three A^^ > and are all marginally irrelevant. In 
this case, the phase is C0S2[lcr, 2a] with two gapless modes 
in the spin sector. SU(2) spin invariance fixes the Luttinger 
parameters in the spin sector, gi^ = g2a = 1- After pinning 
the ^p-^- and ipp^ , the scaling dimensions for the observables 
are 
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pure imaginary , (64) 



Thus we have spin and bond energy correlations oscillating 
with period 4 and decaying with power law 1/x. 



IV. SPIN-GAPPED PHASES IN ORBITAL FIELD 

Besides the spin-gapless phase, C0S2, there are other cases 
in which the spin sector is partially or fully gapped. Below we 
discuss each case in detail and summarize the main properties 
in Table HD 



A. A?i < 0, AJa > 0, A?2 > 

In this case, only Afj^ is relevant and flows to strong cou- 
pling. We pin Oiij such that cos (2-\/20ict) = 1 and the phase 
is COS 1 [2a]. We have S^,,^^ ~ sin(^/202.) and S^'^g ^ 
cos(-\/26'2cr), SO both show 1/x power law correlations. 



B. A?i > 0, A^2 < 0, A?2 > 

In this case, we pin 62^ such that cos (2\/2^2cr) — 1- The 
phase is COSl [Icr] and is qualitatively similar to the previous 
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FIG. 3: Top; Period 4 translational symmetry breaking when 
-^11, -^22 > 0, A12 < 0, drawn in the ID chain picture. The bond 
energy pattern is given by Eqs. J70t-ll7U and the chirality pattern by 
Eq. J72t . Thicker lines represent stronger bond; "+" and "-" sym- 
bols of varying boldness schematize the scalar chirality associated 
with sites (or equivalently with site-centered loops); and arrows on 
the links show the bond currents. Bottom: The same pattern in the 
two-leg triangular ladder drawing. 



C. A?i,Af2>0,A?2<0 

In this case, Afj^ and A22 are marginally irrelevant while 
Aj2 is marginally relevant and flows to strong coupling. To 
minimize the energy associated with AJ2' cf. Eq. ( [38] l, we pin 

9a+ and tp„_ to satisfy. 
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To characterize the resulting COSO fully gapped phase, we 



note that Ej and B^^l„„ gain expectation values, 
late the first- and second-neighbor bond energies. 



We calcu- 



' ^ "-^stagg 

~ A cos \—x+ — +a\ +C COS (ttx) , (70) 



5B^'^\x) - e^t^e*^ (^P'S| + 4^'S|) + 



A cos I —X + — + a 



(71) 



where A, C, and A' are some non-universal real numbers, 
while a is fixed to one of the values {±j, ± ^}. We see 
that this phase has translation symmetry breaking with period 
4 as illustrated in Fig. |3] The four independent values of a 
correspond to four translations of the bond pattern along x. 

To further characterize the state, we also calculate the scalar 
chirality, 

Xix) ^ A" cos (-X -- + a] +C" cos (nx) , (72) 

where A" and C" are some non-universal real amplitudes, 
while a is the same as in Eqs. JTOli-JTni. The period-4 pat- 
tern induced in the chirality is also shown in Fig. [3] and is 
consistent with the spontaneous period-4 bond order on top of 
the staggered chirality background present from the outset. 



D. A?i,A?2 <0,A?2 >0 

In this case, AJj and A22 are marginally relevant and flow 
to strong coupling while AJ2 is marginally irrelevant. To min- 



-H - -I- 



-I- 



-I- - + 



''o^'^^''d^/d^'''Q^''d^''0\ 



FIG. 4: Top: Static period-2 VBS when A?i, A^2 < 0, A?2 > 0, 
drawn in the ID chain picture. Note that the background staggered 
chirality is present from the outset due to the orbital field. Bottom: 
The same pattern in the two-leg triangular ladder drawing. 



imize the relevant interactions, we pin 

cos (2^/26*1^) = cos (2726*2^) = 1 



(73) 



This is a different COSO fully gapped phase where only 
_Bg"~g gain expectation values. The nearest-neighbor bond 
energy is 



SB'^^\x) 






stagg 



C'cos(7ra;) . 



(74) 



The physical picture of this phase is shown in Fig. |4] 



E. \'i2 < and either AJj < or A^2 < 

Here, we do not know how to minimize the relevant interac- 
tions due to the competition of the pinning conditions in Vx, 
Eq. (l38T l. However, we expect that, depending which terms 
grow faster under the RG and win, the final outcome reduces 
to one of the phases discussed above. 



V. DISCUSSION 

In this paper, we considered the effects of orbital field on 
the half-filled electronic two-leg triangular ladder. In weak 
coupling, the Umklapp H" [Eqs. ( fT3] l and ( |32] |1 always makes 
the system Mott-insulating, and we described in detail possi- 
ble phases. 

We would like to conclude by indicating a connection with 
the Spin Bose-Metal (SBM) theory in Ref. |3 and discussing 



effects of the orbital field on the SBM. It turns out that our 
present electronic results translate readily to this case. The 
SBM can be viewed as an intermediate coupling Cl[p— ]S2 
phase and is obtained in the absence of the field by gapping 
out the overall charge mode using an eight-fermion Umklapp 
term, whose bosonized form is^ 



H^ = 2w8COs(46'p+) 



(75) 



Reference argued that wg > is appropriate for the elec- 
tronic case that corresponds to a spin- 1/2 system with ring 
exchanges on the zigzag ladder. This gives pinning condition 
for the overall charge mode. 
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P+ 



(mod 27r) 



(76) 



Note that this pinning condition is compatible with the pinning 
Eq. (|40] | due to the new four-fermion Umklapp "H" arising in 
the presence of the orbital field, so the two Umklapps lead to 
similar Mott insulators. 

We can consider situations where the main driving force to 
produce Mott insulator is the eight-fermion Umklapp while 
the orbital field is a small perturbation onto the SBM phase. 
Formulated entirely in the spin language, the underlying elec- 
tronic orbital fields give rise to new terms in the Hamiltonian 
of a form Si ■ [S2 x 63] on each triangle circled in the same 
direction.—— In the ID chain language, this becomes a stag- 
gered spin chirality term {—lYS{x — 1) • [S{x) x S{x + 1)]. 
Starting from the SBM theory in the absence of the field, this 



gives a new residual interaction of the same form as "H" (sim- 
ilar to Xtt in Ref. |2|). In principle, this T-L^ can be irrelevant in 
the SBM phase if the one Luttinger parameter gp_ in the SBM 
theory^, is less than 1/2, and in this case the orbital effects will 
renormalize down on long length scales. On the other hand, if 
this terms is relevant and pins tpp_ , then the resulting phases 
are precisely as already considered in the electronic language. 
In this simple-minded approach, all the phases we discussed 
in the paper are proximate to the SBM phase. It would be 
interesting to explore spin models realizing the SBM in the 
presence of such additional chirality termsi^i^ 

The presented physics appears to be rather special to the 2- 
leg ladder case, but is quite interesting in the context of such 
models. Perhaps the most intriguing finding is the C0S2 phase 
with two gapless spin modes. Note that the relevant chirality 
interaction involves both chains and the system is far from the 
regime of decoupled chains. Our characterization of this state 
comes from the formal bosonization treatment, but it would 
be interesting to develop a simpler intuitive picture. 
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